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Abstract 

Central heavy-ion collisions may induce sizeable fluctuations of the 
topological charge. This effect is expected to distort the dispersion 
relation for the hadron masses. We construct a general setup for a 
compact description of this phenomenon in the framework of bottom- 
up holographic approach to QCD. A couple of soft wall holographic 
models are proposed for the vector mesons. The states having different 
circular polarizations are shown to have different effective mass. The 
requirement of stability imposes strong constraints on the possible 
choice of models. 


1 Introduction 

Recently the study of light mesons by means of a bottom-up holographic 
approach has attracted a lot of attention(see, e.g., the short surveys [Tj). 
Usually the spectra of resonances and the related physics are analysed at 
vacuum conditions. However, in view of current experiments with relativis¬ 
tic heavy ion collisions at RHIC, GSI, and CERN, it is useful to include in 
the holographic models some non-trivial external conditions. One of theo¬ 
retical methods consists in incorporation of the Chern-Simons (CS) term in 
a holographic action. This permits to address holographically such prob¬ 
lems as the magnetic susceptibility of the quark condensate [2j, the chiral 
magnetic effect mm, some subtle questions in the behavior of the correla¬ 
tion functions M, and derivation of the 0(p 6 ) Chiral Perturbation Theory 
low-energy constants [7]. 

ft has been demonstrated recently [8] that a suitable CS term may lead to 
interesting meson phenomenology. This CS term is motivated by the possi¬ 
bility of local parity breaking taking place in baryonic matter; a phenomenon 
that could be triggered by large topological fluctuations taking place in cen¬ 
tral heavy ion collisions. For light quarks in a quasi-equilibrium situation the 
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creation of a topological charge translates immediately into the generation 
of a finite chiral chemical potential 

Inspired by these observations, in the present paper we analyse the impact 
of the CS term in the background considered in Ref. [8] on the mass spectrum 
of the vector mesons in the Soft Wall (SW) holographic model |9j. This will 
give a generalization of results of Ref. [8] to higher radial excitations in the 
large- N c limit of QCD. Our choice of the model is motivated by its nice 
property of possessing the linear Regge like spectrum which is expected in 
the first approximation both experimentally [10| and in the string like models 
of hadrons. To make the consideration clear we will restrict ourselves by the 
simplest version of the SW model. 

The paper is organized as follows. The general holographic setup is pre¬ 
sented in Section 2. The analysis of impact of the axial chemical potential on 
the mass spectrum is given in Section 3. A couple of exactly solvable models 
are constructed in Section 4. A short discussion of our results are contained 
in the concluding Section 5. 


2 The holographic setup 

We consider the gauge SW model [9] with 5D Abelian fields L and R dual (on 
the AdS 5 boundary) to the sources of the left and right 4D vector currents. 
A parity-odd 5D CS term is added to the action, 


S — S'freefL] + Sf ree [.R] + S'cs[-b] — Scs[A], (1) 

Sfree[B] = J d A xdze^B MN B MN , B — L, R, (2) 

Scs[B] = -kJ d 4 xdze MNARC B M B NA B RC . (3) 

The 5D space is the AdSs one with the metric 

r2 

ds 2 = —(dx^dx^ — dz 2 ), /x = 0,1,2,3, (4) 


where R is the radius of the AdS 5 space and z > 0 represents the holo¬ 
graphic coordinate. The dilaton background e v is not yet fixed for general¬ 
ity. As usual, the constants g§ and k are fixed by matching to the ultraviolet 
asymptotics of the two-point vector correlator m and to the axial anomaly, 
respectively 


9l 12t r 2 N c 

R N c ’ 24vr 2 ‘ 


( 5 ) 
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In terms of the vector, V — (L + R)/ 2, and axial-vector, A = (R — L)/ 2, 
fields the free and CS parts of the action can be rewritten as (Lorentz indices 
are lowered) 



Sqs — —•k J cRxdze MNABC A m ( V NA V B c + A na A bc ) ■ (7) 

If one wishes to provide conservation of the 4D vector current, the Bardeen 
surface counterterm must be added, 

S B = 2kJ d 4 xe^ xp A lj y i ,V Xp . ( 8 ) 

Then one obtains the standard result for the covariant anomaly m , 

= 0, d p J A = :ikV lw V] lu + kA^A^, (9) 

where F ^ = \e pvXp F Xp . 

In contrast to the finite density effects, the thermal corrections to the 
meson masses appear in the next-to-leading order in the large-lV c counting 
as they emerge due to the pion loops. On the other hand, the holographic 
approach is inherently large-iV c one. For this reason we will not include the 
finite temperature effects into our considerations. 


3 Embedding the axial chemical potential 

In Ref. [8] it was assumed that the axial chemical potential /i 5 arises from 
a time-dependent but spatially homogeneous background of a pseudoscalar 
field a(t) such that /i 5 = a(t). In the 5D setup, we will treat A M as a 
background axial-vector field and relate a(t) to the z-component of Am- 
Namely, we assume for the vacuum expectation value of the vector field that 

(A m ) = (A z ) = ix 5 x 0 f{z ), (10) 

where the shape function f(z) will be specified later. The expressions ([9]) 
and (ITU lead to the following form for the vector part of the action (ITT) . 

S = ±f d A xdz (-f °-Vm N + ^5 fe 05ABC V A d B V c ) . (11) 

Here ^ = 1 (see (jSJ)). 
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( 12 ) 


In the axial gauge V z = 0, the equation of motion reads 


d z 



-—d 2 V^- 2n 5 fe mik diV k 
z p 


0. 


The small latin indices denote the usual space coordinates, m,i,k = 1,2,3. 
Making the 4D Fourier transform, V^(x,z) = f d 4 pe ipx V M (p, z) and assuming 
the standard plane wave ansatz, V^(p, z) = e^v^z), we arrive at the equation 
for the particle-like excitations, 

a, (f a,v) + f P \ 

The physical spectrum is given by the eigenvalues p 2 n = rn 2 of normalizable 
solutions of Eq. Et However, the last term induces the mixing between 
different polarizations. We must find a basis diagonalizing the Eq. (Eli. 

For convenience, let us introduce the notation for the differential operator 


£-fj, T i2£ m i k pi£ k v 0. (^3) 


F = d z \ —d z + — p . 

z z 


(14) 


The space-like part of Eq. (ITS]) can be rewritten in the vector form 

[Pe + i2p 5 fp x s^j v — 0. (15) 

The equation (1T5D is diagonalized with the help of the projectors, 

(16) 
P~ 

1 " 


P = 
lk 2 




ik 


p 

PiPk 

p 2 


± 


IpI 


^iknPn 


( 17 ) 


The projectors (1T7T) on the ’’circular” polarizations have the following evident 
properties, 

P ± P^ = 0, p ± p ± = p ± , trP ± = l, P ± p = 0, P + +P-=P ± .( 18) 


Now we change the basis for the space-like polarizations from e = {£\,£ 2 i£p 
to (£||, £_, £ + ). The equation (1131) takes the form 

Fv 11 = 0, (19) 

(P±2p 5 f\p\yp = 0. ( 20 ) 

Thus, the longitudinal and circular polarizations will have different masses, 
with the latter being dependent on the momentum. In other words, instead 
of a ’’peak” corresponding to a given meson we will see three ’’peaks” with 
the splitting depending on the value of three-dimensional momentum. A 
similar phenomenon in a different situation was obtained in Refs. 0Q3]- 
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4 Solvable models 


We have not yet specified the dilaton background <p(z) and the shape func¬ 
tion f(z) for the axial chemical potential. Various holographic models for 
the ’’splitting phenomenon” can be obtained by fixing these functions. The 
simplest solvable SW model resulting in a Regge-like spectrum is given by 
the background [9]. 

(p = -A V. (21) 

Let us accept this background and consider the longitudinal polarization. 
Following Ref. [9], we use the substitution 

v(z) = e x2z2 ' 2 ^(z), (22) 

to convert the Eq. (1191) into the Schrodinger form 

- d 2 M + (aV + = mfyl [. (23) 


Introducing the dimensionless variable y = \z, Eq. f!23p for the discrete mass 
spectrum m 2 = p ^ of the longitudinal (and time like) polarization transforms 
into 

( o \ TTl 2 

y2 + i’n = (24) 

The normalized solutions are (n — 0,1, 2,...) 


V 1 ! = 


2rd 


(1 + n) 




(25) 


where L\ denote the associated Laguerre polynomials. The mass spectrum 
is given by the corresponding eigenvalues, 


m l,l = 4A2 (™ + !)> 


(26) 


where the parameter A controls the slope of the radial Regge trajectory. 

To obtain the spectrum of circular polarizations we need to specify the 
function f(z). The choice of f(z) must comply with the requirement of 
correct UV/IR behavior of the action with respect to conformal symmetry. 
This constraint leaves, however, much freedom. We will be interested in the 
exactly solvable cases. A simple possibility of this sort is given by the ansatz 


/ 


b e~ x2z2 
2 z 


(27) 
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Here b is a dimensionless constant. After the substitution (1221) the Schrodinger 
equation (124)) acquires an additional contribution, 

-d 2 y^t+ (y 2 + ^ ± ^-\P\')^n ^ ^T^n- ( 28 ) 

The eigenfunctions remain the same as in (1251) but the mass spectrum is 
shifted, 

m 2 n± = 4A 2 {n + 1) ± bp 5 \p\. (29) 

Thus, the massive vector fields split into three polarizations with masses 
m n - < m n) || < m n ^ + . The mass splitting is linearly dependent on the value 
of the spatial momentum p. The formula (1291) can be considered as a gener¬ 
alization of the result of Ref. [8] to the radially excited spectrum. 

The matching of (125|) with the corresponding expression in Ref. [8], 

m v,± — m v =•= C|p|, (30) 

allows to estimate the constant b. Within the effective model studied in 
Ref. [8], ( = N c g 2 fi 5 /87i 2 , where g p is related with the mass of vector particle, 
m 2 = 2g 2 f 2 ~ rri^. Substituting the phenomenological values for m p and for 
the weak pion decay constant f n , one obtains the estimate ( ~ 1.5/i 5 . We 
remark that this estimate can be made without phenomenological values if 
one uses the formula for rri p from the QCD sum rules in the large-iV c limit 
(see, e.g., nap, rn 2 = 247 r 2 f 2 /N c that yields g p = I2n 2 /N c and leads directly 
to ( = |/i 5 . Comparing (129|) and (1301) we arrive at the value b — |. 

Using the results of Ref. [15], the background (1211) can be generalized to 

cp = — X 2 z 2 log U 2 (c, 0; X 2 z 2 ), (31) 


where U is the Tricomi hypergeometric function. Then the dimensionless 
parameter c in (l3Tj) will control the shift of the radial Regge trajectories: 
n + 1 in (1261) and (1291) will be replaced by n + 1 + c. 

The model was designed such that the relation m 2 + — _ = const holds 

for any n. Other possibilities are of course possible. For instance, one can 
construct a solvable model with the relation rn 2 n+ /m 2 n _ = const. This is 
achieved via replacing (1271) by 



Introducing the variable 


V = 


i±^rlp| 


A z, 


(32) 


(33) 
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the analogue of Eq. (1251) looks as follows 


dfyn + ( y z + ) C = 




L± 


m, 


n,± I ± 

t'n, 


AVl±^|p| 


resulting in the mass spectrum 


bp 5, 


m n,± = 4A \ll±—\p\ (ra + 1 ). 


(34) 


(35) 


The eigenfunctions (1251) are now momentum- and / 15 -dependent because the 
variable y (1551) depends on \p\ and / 15 . The matching to Eq. (130|) for n — 0 
for small \p\ gives b = ^A 2 . In this scenario, the contribution of /15 to the 
masses grows with n. 

Since the effective masses of polarized modes depend on the spatial mo¬ 
mentum in the models under consideration, a certain care must be exercised 
to provide stability and, as a consequence, to avoid the superluminal propa¬ 
gation. From (129]) we have for the energy 

Po,n = P 2 + 4A 2 (nT 1) ±bp 5 \p\. (36) 


In this respect, the £_ polarization is dangerous. Imposing the condition 
Po n — 0 we arrive at 

I « \2 

/i§<-y-(n+l). (37) 

The strongest limitation comes from the ground state n — 0, 

4A , . 

hs < y. (38) 

Consider now the group velocity v = ^. It cannot exceed the speed of 
light in the vacuum, v < 1. For (j55j) this yields 


_ 2 jgj ± bfi 5 _ 

2a Jp 2 + 4A 2 (n + 1) ± bfi 5 \p\ 


< 1 . 


(39) 


The constraint (1521) leads to the same limitations (1571) and (1551) . 

The analogue of constraint (1591) for the second model considered above 
takes the form 


IpI ± 


6/^5 ( re + l ) 


\ 2 y/ l±b/i5|p|/A 4 


< 1 . 


p 2 + 4A 2 \/l ± 6/x 5 |p|/A 4 (n + 1) 


(40) 
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It is seen that the condition (HUD cannot be fulfilled at any n — the larger 
is n the stronger is limitation on /i 5 . Since n labels the Kaluza-Klein modes 
of a five-dimensional field, all infinite number of resonances must be present 
simultaneously in agreement with the limit N c —> oo. Thus, the second 
solvable model should be disregarded. 

5 Concluding discussions 

According to some estimates of Ref. [8], the natural value of /i^ in the central 
heavy ion collisions is in the few hundreds of MeV. Such values are far below 
the limitation (1381) . 

/r 5 < ^ ^ « 1400 MeV, (41) 

0 o 

where the value A ~ 530 MeV is taken from the typical slopes of the vec¬ 
tor radial trajectories. Thus, the first model passes the stability criterium. 
We have shown that this criterium imposes a strong constraint on the pos¬ 
sible choice of models. In particular, it seems to falsify the second model 
constructed in the previous Section. 

In principle, the effect of mass splitting between different polarizations can 
be quite large. The masses of £ + and £_ circular polarizations of neighboring 
resonances can even overlap. 

We expect that the contribution of the axial chemical potential to the 
hadron masses is approximately constant, i.e. it does not depend on the 
radial number n. In this regard, the first constructed model looks physical. 

The splitting of effective masses of two polarizations signifies the parity 
breaking in the medium. In the framework of Ref. [8], this effect is related 
with the violation of 4D Lorentz invariance due to the time-dependent back¬ 
ground which is present in the Lagrangian from the very beginning. In the 
presented holographic description, we have essentially the same situation. 

The experimental implications of the presence of the axial chemical po¬ 
tential are discussed in detail in Ref. [TU j. 

Finally we mention some possible future directions. First of all, it would 
be interesting to study the impact of the axial chemical potential on the 
masses of vector mesons in the top-down holographic approach. Second, it 
makes sense to include a non-zero isospin chemical potential /ij into our 
consideration. In the low-energy QCD, jU/ contributes to the masses of 
charged pions, m T ± ~ m n o ± /ij, and triggers the pion condensation at val¬ 
ues fii > m^o D3- This result is reproduced by the hard-wall bottom-up 
holographic model [18] and by the Sakai-Sugimoto top-down holographic 
model p(9]. In addition, the latter analysis shows that at jii > 1.7m p , the 


lowest vector meson also condenses. It is thus interesting to verify this effect 
in the bottom-up approach and then to study the interplay of [ij and the 
axial chemical potential. 
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